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Abstract 

The twist and writhe numbers and magnetic energy of an orthogonally perturbed vortex 
filaments are obtained from the computation of the magnetic helicity of geodesic and ab- 
normal magnetohydrodynamical (MHD) vortex filament solutions. Twist is computed from 
a formula recently derived by Berger and Prior [J. Phys. A 39 (2006) 8321] and finally 
writhe is computed from the theorem that the helicity is proportional to the sum of twist 
and writhe. The writhe number is proportional to the total torsion and to integrals of the 
vector potential. The magnetic energy is computed in terms of the integral of the torsion 
squared which allows us to place a bound to energy in in the case of helical filaments due 
to a theorem by Fukumoto [J. Phys. Soc. Japan, (1987)] in fluid vortex filaments. It is 
also shown that filament torsion coincides with the magnetic twist in the case under con- 
sideration, where a small orthogonal magnetic field exist along the thin filament. A new 
Aharonov-Bohm (AB) phase term is obtained in the writhe number expression which is not 
present in the Moffatt-Ricca (Proc Roy Soc London A, 1992) computation. PACS numbers: 

02.40 
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I Introduction 



The topology and geometry of fluid vortex filaments [1] have immensely contributed to our 
comprehension of the [2] dynamics and kinematics of magnetic filamentary structures with 
applications in solar and plasma physics [3, 4]. In particular twisted filamentary structure have 
helped solar physicists to better understand the mechanism which allow the highly twisted 
coronal magnetic flux tube to emerge from the solar surface and produce those beautiful solar 
flares and loops. Recently Berger and Prior [5] have presented a detailed discussion of twist 
and writhe of open and closed curves presenting a formula for the twist of magnetic filaments in 
terms of parallel electric currents. In the present paper we apply Berger-Prior formula and the 
theorem of the sum of writhe and twist numbers investigated previously by Moffat and Ricca 
[6] and Berger and Field [7] to obtain an expression for the writhe number for the magnetically 
perturbed vortex filamentary twisted structure. To obtain the writhe number we make use 
of the helicity local expression for the magnetic field proportional to the magnetic field itself 
and decompose the magnetic vector field along the Serret-Frenet frame in 3D dimensions. 
This allows us to solve scalar MHD equations to obtain constraints on the global helicity 
expression which will allow us to compute the writhe number. The knowledge of tilt , twist 
and writhe of solar filaments for example has recently helped solar physicists [8] to work out 
data obtained from the vector magnetograms placed in solar satellites. This is already an 
strong motivation to go on investigating topological properties of these filamentary twisted 
magnetic structures. Earlier the Yokkoh solar mission has shown that the sigmoids which 
are nonplanar solar filaments are obtained due to the action of electric currents along these 
filaments which was used recently [9] as motivation to investigate current-carrying torsioned 
twisted magnetic curves Throughout the paper we use the notation of a previously paper on 
vortex filaments in MHD. Mathematical notation is used based on the book by C. Rogers and 
W. Schief [10] on the geometry of solitons. The paper is organized as follows: In section 2 
we decompose MHD equations on a Frenet frame along the twisted thin filament and solve 
the scalar equations obtained. In section 3 we compute the twist, writhe and energy of the 
magnetic filament. In section 4 we present the conclusions. 
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II MHD scalar equations for twisted filamentary struc- 
tures 



Let us now start by considering the MHD field equations 

V.S = (II.l) 

VxB = aB (11.2) 

where a is the magnetic twist and the magnetic field B along the filament is defined by the 
expression B = B s t+B n t and B n is the magnetic field perturbation orthogonal to the filament 
all along its extension, and B s is the component along the arc length s of the filament. The 
vectors t and n along with binormal vector b together form the Frenet frame which obeys the 
Frenet-Serret equations 

? = nn (II.3) 

ft = -Kt + Tb (II.4) 

b' = -rn (II.5) 

the dash represents the ordinary derivation with respect to coordinate s, and «(s, £) is the 
curvature of the curve where k = Rr 1 . Here r represents the Frenet torsion. We follow the 
assumption that the Frenet frame may depend on other degrees of freedom such as that the 
gradient operator becomes 

->d d -> d 
V = t Fs +n ml + b 8i (H6) 

The other equations for the other legs of the Frenet frame are 

d - 

—t = e ns n+[Vt b + r]b (11.7) 
on 

d 

—n = -9 ns t - (divb)b (11.8) 
b=-[n b + r]t- (divb)n (11.9) 



dn 
~db 



^-t = e bs b-[n n + r]n (11.10) 



—n = \n n + r]t- k + (divn)b (II. 11) 

do 



d -> 

—b = -6 b8 t - [k + (divn)]n (11.12) 
do 

Substitution of these equations into the magnetic helicity equation reads 

VxB = i[B n {divb) - B S {Q S + r)] + nrB n + b[nB s - B n (l + 9 ns )} (11.13) 
while its RHS is 

aB = a[B n n + B s t\ (11.14) 
Comparing both sides component by component one obtains the three scalar equations 

rB n = aB n (11.15) 

This result seemed to be noticed earlier by Parker [11] which used to call the parameter a 
torsion. Throughout this derivation we consider that the parameter a is constant, which is 
not very usual in other works on magnetic helicity. The other scalar equations are 

B„ K 



B s (1 + ns ) 

which is also usual when one deals with magnetic flux tubes [12] and finally 



(11.16) 



uB s = B n divb - B S (Q S + r) (11.17) 
The equation V.B = becomes 

d s B s + [6 bs + divb]B s = (11.18) 

A simple solution of this equation can be obtained if one considers that divb does not vary 
appreciably on a short distance ds along the filament. The solution is 

B s = B [1- divb J nds] (11.19) 

where B is an integration constant. Substitution into equation (11.17) assuming that the flow 
is geodesic along the filament and abnormality relation f2 s = along with 6f, s = 9 ns = Oone 
obtains 

B n = BoKdivb)- 1 - J Kds] (11.20) 

which allows us to say that we obtain a solution for the magnetic filament in terms of the 
vortex filament invariant / nds. In the next section we shall compute the magnetic energy in 
terms of the magnetic vector potential A, which allows us to compute the helicity integral and 
the writhe number in terms of the AB phase. 



Ill Magnetic energy, Twist and Writhe 



To be able to compute the magnetic energy , twist and writhe we now compute the vector 
potential assuming that it obeys the Coulomb gauge V.A = and the definition B = VxA 
The Coulomb gauge becomes 

d s A s + divnA n = (111.21) 
together with the equations for the definition of B, taking A = A s t + A n n + A b b yields 

4 1 = -- (HI.22) 

A b K 

which was obtained from the constraint B s = and 

B n = -rA n (111.23) 

B s = - [2tA s + divbA n + (k + divn)A b ] (111.24) 
Algebraic manipulation of the above equations allow us to compute the magnetic energy 

E B = — [ B 2 dV (111.25) 



8tt 

which yields 

Eb= ~LI (jB * 2 + Bn ^ dv (IIL26) 

or in terms in terms of the vector potential component A n yields 

E B = [ J r 2 ds + 1 j divbds] (111.27) 

which upon the assumption divb « r 2 reduces to 

E B = ^-^[[ T 2 ds} (111.28) 
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In case of helical filaments one knows that r = c k where c is constant and the expression 
(111.27) reads 

E B = C -^^[J n 2 ds\ (111.29) 
Due to a theorem by Fukumoto [13] 



J K 2 ds < 16tt 2 L (111.30) 



which clearly places a bound on the magnetic energy E B . Here L denotes the length of the 
filament. Algebraic manipulation of the same equations yields a relation between the Frenet 
curvature and torsion in terms of divn as 

.±V2r = -^ (111.31) 

In the helical case we obtain the following equation 

r-, divn , TTT 
k[1 ±V2] = — (111.32) 

Now let us compute the global helicity 

H = J A.BdV (111.33) 

which is equivalent to 

H = f [A S .B S + A n .B n \dV (111.34) 
substitution of the relations above yields 

H=- J \2rA n + divb]A n ds (111.35) 

Due to a theorem which states 

H= ^-[Tw + Wr] (111.36) 

where $ is the magnetic flux / B.dS. Before computing the writhe number we need to compute 
the twist of the magnetic filament. But this becomes quite simple thanks to a formula recently 
derived by Berger and Prior [5] which is given by 

d m J\\ aBc 



—Tw — — — - = a (111.37) 

ds \B\ B s 



where J\\ is the electric current along the twisted filament. This result was obtained thanks 
to the assumption that B n «< B s since B n is nothing but a simple perturbation of B s . 
Integration of the last expression yields 

Tw = J rds (111.38) 

Substitution of the twist and magnetic helicity H into expression (111.36) yields finally the 
expression for the writhe number 

Air A 2 r 2tt -> r 

Wr = -[(1 + -^-) J rds + ^(divb) J A n ds] (111.39) 

Note that the last term in this expression is proportional to the AB phase / A n ds. 
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IV Conclusions 



In conclusion, a new expression for the writhe of magnetic vortex twisted filament has been 
obtained presenting a new term representing a Berrys phase which was not present in the 
previous calculation of Moffatt and Ricca. A bound in the energy has been obtained thanks 
to a theorem by Fukumoto on the bound of the integral invariant of the square of the Frenet 
curvature value in the helical case. Future applications in plasma physics or in DNA [12] 
be appear elsewhere. Expressions for the writhe number for the elastodynamics have been 
obtained by Klapper and Tabor [13]. Other integral invariants in vortex fluid dynamics have 
been recently also computed by Maggioni and Ricca [14] which can be generalised to MHD. 
Yet other interesting examples of the topological bounds to energy [2] have been provided by 
Khesin [15]. 
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